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Abstract. We present the first high-resolution (5'.' 6 x 5'.' 1) images of the emission of silicon monoxide (SiO) in the nucleus of 
the nearby spiral IC 342 obtained with the IRAM Plateau de Bure Interferometer (PdBI). Using a two-field mosaic, we have 
mapped simultaneously the emission of the SiO(v = 0, J = 2- 1) and H I3 CO + (7 = 1 -0) lines in a region of ~0.9kpcx 1.3 kpc 
(RAxDec)-size centered around the nucleus of IC 342. The bulk of the emission in the two lines comes from a ~ 290 pc spiral 
arm located North and a central component that forms the southern ridge of a r ~ 80 pc nuclear ring which has been identified 
in other interferometer maps of the galaxy. We detect continuum emission at 86.8 GHz in a ~ 80 - 180 pc central source. The 
continuum emission, dominated by thermal free-free bremsstrahlung, is mostly anticorrelated with the observed distribution of 
SiO clouds. The SiO-to-H 13 CO + intensity ratio is seen to increase by an order of magnitude from the nuclear ring (~ 0.3) to the 
spiral arm (~ 3.3). Furthermore the gas kinematics show significant differences between SiO and H 13 CO + over the spiral arm 
where the linewidths of SiO are a factor of 2 larger than that of H 13 CO + . The average abundance of SiO in the inner r ~ 320 pc 
of IC 342 is X(SiO)> 2 x 10~ 10 . This evidences that shock chemistry is at work in the inner molecular gas reservoir of IC 342. 
To shed light on the nature of shocks in IC 342, we have compared the emission of SiO with another tracer of molecular shocks: 
the emission of methanol (CH3OH). We find that the significant difference of the abundance of SiO measured between the spiral 
arm (X(SiO)~a few 10~ 9 ) and the nuclear ring (X(SiO)~10~ 10 ) is not echoed by a comparable variation in the SiO-to-CH^OH 
intensity ratio. This implies that the typical shock velocities should be similar in the two regions. In contrast, the fraction of 
shocked molecular gas should be ~ 5 - 7 times larger in the spiral arm (up to ~10% of the available molecular gas mass over 
the arm region) compared to the nuclear ring. In the light of these results, we revise the validity of the various scenarios which 
have been proposed to explain the onset of shock chemistry in galaxies and study their applicability to the nucleus of IC 342. 
We conclude that the large-scale shocks revealed by the SiO map of IC 342 are mostly unrelated with star formation and arise 
instead in a pre-starburst phase. Shocks are driven by cloud-cloud collisions along the potential well of the IC 342 bar. The 
general implications for the current understanding of galaxy evolution are discussed. 

Key words, galaxies: individual: IC 342 - galaxies: starburst - galaxies: nuclei - ISM: molecules - molecular processes - 
radio lines: galaxies 



1. Introduction 

There is mounting evidence that the properties of molecular 
gas in starbursts (SBs) and Active Galactic Nuclei (AGNs) dif- 
fer from that of quiescent star forming galaxies (e.g., Genzel et 
al. 1998). The spectacular energies injected into the gas reser- 
voirs of active galaxies can create a particularly harsh envi- 
ronment for the neutral ISM. Although thus far restricted to 
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a handful of objects, multiline millimeter wave studies have 
allowed to study the onset of large-scale shocks, the prop- 
agation of chemistry of Photon Dominated Regions (PDRs) 
or the prevalence of X-ray Dominated Regions (XDRs) in 
the molecular disks of SBs and AGNs (e.g., Mauersberger & 
Henkel 1993; Tacconi et al. 1994; Garci'a-Burillo et al. 2000, 
2001b, 2002; Martin et al. 2003, 2005; Usero et al. 2004; 
Fuente et al. 2005a; Meier & Turner 2005). 

The first SiO(v = 0, 7 = 2-1) maps made with the Plateau 
de Bure Interferometer (PdBI) in the nuclei of the prototypi- 
cal starbursts NGC 253 and M 82 have revealed the existence 
of large-scale molecular shocks in galaxy disks (Garcia-Burillo 
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et al. 2000, 2001b). Different scenarios have been proposed to 
account for the emission of SiO in our own Galaxy and in the 
nuclei of external galaxies. On small ~pc-scales, studies in the 
Galaxy disk show that the enhancement of SiO in gas phase 
can be produced in the bipolar outflows of young stellar objects 
(YSOs), due to the sputtering of dust grains by shocks (Martin- 
Pintado et al. 1992; Schilke et al. 1997; Cesaroni et al. 1999). 
On larger scales, Martin-Pintado et al. (1997) reported the de- 
tection of a SiO ~150 pc circumnuclear disk (CND) in the 
Galactic Center region. In this CND high fractional abundances 
of SiO are found in molecular clouds which are not actively 
forming stars, but where bar models for our Galaxy predict a 
high likelihood for cloud collisions (Huttemeister et al. 1998). 
In M 82, virtually all of the SiO emission traces the disk-halo 
interface where episodes of mass injection are building up the 
gaseous halo (Garcia-Burillo et al. 2001b). Garcia-Burillo et 
al. (2000) have discussed the role of bar resonances at induc- 
ing shocks in the -600 pc CND of NGC 253. However, the 
high-inclination of NGC 253 and the limited perspective of 
the Galactic Center region make the determination of bar reso- 
nance positions mostly dependent on kinematical models. 

In this paper we study at high-resolution (~5") the SiO(v = 
0, 7 = 2-1) emission in the inner r ~ 30" of IC 342, using 
the IRAM Plateau de Bure Interferometer (PdBI). IC 342 is a 
nearby (D^3.3 Mpc, i.e., 1"^16 pc; Saha, Claver & Hoessel 
2002) weakly barred spiral galaxy which hosts a moderate star- 
burst episode in the central r ~ 80 pc nuclear region (Boker et 
al. 1997, 1999). Thanks to its nearly face-on orientation and 
close distance, IC 342 is an optimal testbed where the underly- 
ing mechanisms of large-scale molecular shocks can be probed 
using the PdBI. At this distance, the PdBI can trace and spa- 
tially resolve the SiO emission of shocked molecular gas on 
scales of individual GMCs (~ 80 pc) in the disk of IC 342. 
Several works have underlined the similarities between the nu- 
cleus of our Galaxy and that of IC 342 in terms of the measured 
gas mass fractions, stellar masses and star formation (Downes 
et al. 1992). This study can thus help to shed light on the origin 
of molecular shocks in the nucleus of our own Galaxy. 

The distribution of molecular gas in the inner r ~ 320 pc 
of IC 342, revealed by the published interferometer CO maps 
of the galaxy, is reminiscent of the typical response of gas to 
a bar potential (Lo et al. 1984; Ishizuki et al. 1990; Levine et 
al. 1994; Meier & Turner 2001, 2005; Schinnerer et al. 2003). 
Two gas lanes are shifted with respect to the major axis of 
the ~ 9 kpc bar oriented with a position angle ~ 20° (Buta 
& McCall 1999). The gas lanes delineate a two-arm spiral pat- 
tern that ends at a nuclear r ~ 80 pc ring. This pattern would 
correspond to the transition from Xi orbits (outer disk) down to 
X2 orbits (inner disk) of the IC 342 bar, assuming that an Inner 
Lindbland Resonance (ILR) exists near r ~ 80 - 160 pc. 

Previous interferometer maps have probed the dense gas 
content of the nucleus of IC 342 (Ho et al. 1990: NH 3 ; Downes 
et al. 1992: HCN; Nguyen-Q-Rieu et al. 1992: HCO + ; Meier 
& Turner 2005: HNC, HC 3 N, C 2 H, C 34 S, HNCO, CH 3 OH and 
N 2 H + ). Five major GMCs, labelled originally as A-to-E in the 
HCN map of Downes et al. (1992) and later redefined by Meier 
& Turner (2001), are identified in the inner r ~ 320 pc disk. 
Meier & Turner (2005) find remarkable differences in morphol- 



ogy between the various molecular emission maps of IC 342. 
This is interpreted as an evidence of strong chemical differen- 
tiation in the nucleus of IC 342. While some molecules trace 
the Photon Dominated Regions close to the nuclear starburst 
(r ~ 80 - 180 pc), the emission of molecular species such as 
methanol (CH3OH) is considered to be stemming from shocks 
(Meier & Turner 2005). 

In the scenario of shock chemistry there is evidence that 
CH3OH and SiO trace distinctly different velocity regimes in 
shocks. SiO is seen to be associated with more energetic events, 
i.e., those potentially more efficient at processing dust grains 
(Garay et al. 2000). In this paper we take advantage of the com- 
plementarity of SiO and CH3OH as tracers of shock chemistry 
and use the SiO-to-CHsOH ratio in IC 342 to discuss the ori- 
gin of large-scale shocks in this galaxy. With this aim we make 
a quantitative comparison of our results with those obtained in 
other well-known references for shock chemistry in our Galaxy 
and in external galaxies. 

We describe in Sect. 2 the PdBI observations used in this 
paper. Sect. 3 presents the main results issued from the analy- 
sis of the continuum image and the SiO and H 13 CO + (1-0) line 
maps of IC 342. We derive in Sect. 4 the fractional abundances 
of SiO in the disk of the galaxy and interpret their relation with 
the published CH3OH map in Sect. 5. Sect. 6 discusses the 
possible mechanisms driving large-scale shocks in IC 342. In 
Sect. 7 we analyze the potential role of density waves as drivers 
of the large-scale shock chemistry in IC 342. The main conclu- 
sions are summarized in Sect. 8. 



2. Observations 

Observations of IC 342 were carried out with the PdBI from 
July to August 2001. We observed simultaneously the J=2- 
1 line of SiO (86.847 GHz) and the J=l-0 line of H 13 CO + 
(86.754 GHz) using the CD set of configurations. The pri- 
mary beam of the PdBI at 87 GHz is 55". Two positions 
shifted (0",-12") and (0",+12") from the phase center, 
aj2000.o=03 h 46 m 48 s :01, <5 J2 ooo.o=68°05'4670, were observed in 
mosaic mode. We adjusted the spectral correlator to give a 
contiguous bandwidth of 1500 km s -1 . The frequency resolu- 
tion was set to 1.25 MHz (4.3 km s -1 ) during the observations; 
channels were resampled to a velocity resolution of 5 km s~' 
in the final maps. We calibrated visibilities using 0224+671 as 
amplitude and phase reference. The absolute flux scale was de- 
rived on MWC 349, and receiver passband was calibrated on 
3C 454.3 and 3C 345. 

Mosaics were CLEANed using the MAPPING procedure 
of the GILDAS software package, which includes primary 
beam correction. The synthesized clean beam is 576x571 size 
(PA=321°) for the line maps. Images are 300 x 300 pixels in 
extent, with a pixel size of 0735. The rms noise level in 5 km s -1 
wide channels, derived after subtraction of the continuum emis- 
sion, is 1.5 mJy beam -1 at the center of the maps. 

A 3.5 mm continuum map was generated averaging chan- 
nels free of line emission. Uniform weighting was applied to 
the measured visibilities, producing a clean beam of 471 x 378 
(PA = 1 11°). The rms at the center is 0.45 mJy beam -1 . 
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Fig. 1. We overlay the 3.5 mm continuum map obtained with the PdBI 
(contour levels: 0.90 to 4.5 in steps of 0.45 mjy beam -1 ) with the 
Ho-+continuum image (color scale in arbitrary units) obtained by HST 
towards the nucleus of IC 342. Positions of GMCs A, B, C and E 
(as defined by Meier & Turner 2001), are highlighted; the position of 
GMC D, beyond the displayed region of this figure, has been redefined 
in this work according to the local maximum of SiO integrated inten- 
sity (see Fig. 2). A cross marks the position of the galactic center de- 
termined by Schinnerer et al. (2003). Relative positions are referred to 
the phase center ff J2 ooo.o = 03 h 46 m 48 s .01, c5,20oo.o = 68°05'46'.'0. The 
4'.'1 x 3'.'8 beam is drawn at the bottom-right corner. 



3. Results 

3.1 . The 3.5 mm continuum map 

The 3.5 mm continuum emission contours are displayed in 
Fig. 1. The total flux integrated within the inner r ~ 180 pc 
of the galaxy is ~17 mjy. The flux recovered is of ~22 mjy 
when derived from a lower 9" resolution version of the PdBI 
map; the latter is obtained assuming a Gaussian UV taper on 
the visibilities (with a FWHM=70 m). These values are roughly 
in agreement with the previous estimates obtained by Downes 
et al. (1992) and Meier & Turner (2001) at similar frequen- 
cies. The bulk of the continuum emission comes from the inner 
r ~ 80 pc of IC 342 and it is closely linked to the nuclear star 
forming region identified in the Hc+continuum HST image of 
the galaxy (Fig. 1). As illustrated in Fig. 2, the continuum emis- 
sion is mostly anticorrelated with the emission coming from the 
dense molecular gas traced by SiO or H 13 CO + . 

Despite its compactness, the continuum source is spatially 
resolved by the PdBI beam. The morphology of the emission 
follows closely the distribution of star forming complexes in 
the nucleus of IC 342: two 5 Myr-old Hn complexes (Boker 
et al. 1997) which are close to GMCs B and C (notation of 
Meier & Turner 2001), and an older (6-60 Myr) star cluster 
close to the center of the galaxy (Boker et al. 1997, 1999). As 
shown in Fig. 1, two peaks of emission shape the morphology 
of an elongated disk oriented along PA^ 64°. The principal 
peak is related to the western Hn complex, which is identified in 
the HST map to be close to GMC B. The secondary maximum 
is close to the dynamical center where the old star cluster is 
detected. Lower-level emission extends East from the central 



disk towards GMCs A and C. The latter extension is close to 
the eastern Hn complex. 

Accounting for the differences in spatial resolution and sen- 
sitivity, the morphology of the 3.5 mm continuum PdBI map 
agrees with that of previous radiocontinuum maps obtained at 
other wavelengths (Condon et al. 1982: 21 cm; Turner & Ho 
1983: 2 and 6 cm; Ho et al. 1990: 1.3 cm; Downes et al. 1992: 
3.4 mm). At higher frequencies, however, the contribution from 
dust to the thermal emission may not be negligible. This would 
explain the differences between the continuum map of Fig. 1 
and the 1.3 mm continuum image obtained by Meier & Turner 
(2001). 

Based on measurements at various frequencies, previous 
works have concluded that thermal free-free bremsstrahlung 
should dominate the emission budget at 3.5 mm (Downes et 
al. 1992; Turner & Ho 1983; Turner & Hurt 1992). The good 
spatial coincidence between the 3.5 mm and the Ha emissions 
supports the conclusion that the 3.5 mm emission traces the 
location of ongoing star formation in the inner r ~ 80 pc of 
IC 342. 

3.2. The line maps: SiO and H 13 CO + 

3.2.1. Integrated intensity maps and line ratios 

Fig. 2 shows the velocity-integrated intensity maps of SiO(2- 
1) and H 13 CO + (1-0) in the inner r ~ 320 pc of IC 342. The 
spatial resolution of the PdBI map allows to resolve the overall 
emission which, for both tracers, is elongated along the N-S di- 
rection. Though with significant differences between SiO and 
H 13 CO + , the basic morphology of the maps is roughly in agree- 
ment with that emerging from the HCN map of Downes et al. 
(1992). In the case of SiO, we detect strong emission in the 
spiral arm located North (GMCs C and D in Fig. 2). Weaker 
SiO emission delineates the southern ridge of the nuclear ring 
(GMCs A and E in Fig. 2). Finally, SiO emission is not de- 
tected over the southern spiral arm. This result agrees with the 
overall picture derived from other dense gas tracers in IC 342 
which are hardly detected over the southern spiral. Compared 
to SiO, the emission of H I3 CO + is weaker along the northern 
spiral arm, while the peak of emission is found on the southern 
ridge of the nuclear ring. 

As it is shown in Fig. 3, the different distributions of SiO 
and H I3 CO + in IC 342 translate into an order of magnitude 
difference in the SiO-to-H 13 CO + intensity ratio, Ri, which 
goes from ~ 3.3 on the northern spiral arm down to ~0.3 
on the nuclear ring. The average value of Ri inside the im- 
age field of view is ~1.6. If we assume that the emission of 
both lines is optically thin, the value of /?/ provides an esti- 
mate of the beam-averaged fractional abundance of SiO rel- 
ative to H I3 CO + , which is accurate within a factor of 4 (see 
Sect. 4 for discussion). The reference studies of galactic clouds 
(Martin-Pintado et al. 1992; Bachiller & Perez-Gutierrez 1997; 
Fuente et al. 2005b) and external galaxies (Garcia-Burillo et al. 
2000, 2001b; Usero et al. 2004) indicate that a value of /?/ >0.1 
is a strong indication that shock chemistry is at work in molec- 
ular gas. The reported values of R; for IC 342 are similar to 
those found on similar spatial scales in the circumnuclear disk 
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Fig. 2. a(left) We overlay the SiO(2-l) integrated intensity map (contour levels from 90 to 522 in steps of 54 mjy beam -1 km s -1 ; lcr = 
36 mJy beam -1 km s -1 ) with the 3.5 mm continuum map (color scale from 0.90 to 4.5 mjy beam -1 ; lcr = 0.45 mjy beam -1 ), b(right) Same as 
a(left) with H 13 CO + (1-0) (contours) overlaid with SiO(2-l) (color scale), both displayed with the same intensity spacing (used in a(left)). The 
two integrated emission maps were obtained with a 1.2<r clipping applied to the velocity interval [-10, 85] km s -1 . Line contours are screened 
beyond the displayed field-of-view as S/N ratio is lower than 2.5 farther out. Beam-size is represented by a filled ellipse. 



of NGC 253 (-1-3; Garcia-Burillo et al. 2000) and, also, in the 
chimney and the supershell of M 82 (~0.4-3.5; Garcia-Burillo 
et al. 2001b). 

Of particular note, Rj reaches the largest values along the 
northern spiral arm, i.e., in a region of the disk of IC 342 where 
the evidence of active star formation is scarce. Regardless of its 
origin, the order of magnitude variation of/?/ measured over the 
disk of IC 342 reveals that shocks are processing molecular gas 
with a highly changing efficiency (in terms of the total mass 
of grain material processed by shocks relative to the total gas 
mass; see discussion in Sects. 4 and 5). 

3.2.2. Gas kinematics and line profiles 

Fig. 4 shows the velocity-channel maps of SiO(2-l) and 
H 13 CO + (1-0) in the nucleus of IC 342. For both tracers we 
see the expected velocity gradient due to the rotation of the 
disk, which according to the fit of Crosthwaite et al. (2001) 



should be maximal along the kinematic major axis at PA=37°. 
The channel maps show that the gas kinematics are similar for 
SiO and H 13 CO + on the nuclear ring: the velocity centroids 
and linewidths measured in SiO and H 13 CO + are roughly in 
agreement in this region (GMCs A, B, and E). However, the 
gas kinematics show significant differences between the two 
species over the northern spiral arm (GMCs C and D): the emis- 
sion of SiO is detected from v = 16 km s -1 to v = 66 km s -1 
near GMC D, i.e., roughly twice the corresponding velocity in- 
terval for H 13 CO + . 

The different kinematics of the SiO and H 13 CO + lines over 
the northern spiral arm are illustrated by the position-velocity 
plots shown in Fig. 5. SiO lines become significantly wider 
than H 13 CO + lines at the passage of the northern spiral arm 
GMCs (C and D), in contrast with the nuclear ring GMCs 
(B and E) where linewidths are the same within the errors. 
While the velocity centroids derived in the two lines are simi- 
lar in the spiral arm region, the measured linewidths for SiO 
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Fig. 3. SiO(2-l)-to-H 13 CO + (l-0) integrated intensity ratio in the in- 
ner r ~ 320 pc disk of IC 342. The ratio is derived from the two maps 
of Fig. 2, assuming a 2.5cr-clipping on the integrated intensities of 
both lines. The measured ratios (color scale) range from 0.3 to 3.3. 



NIR by Boker et al. (1997) (see Garci'a-Burillo et al. 2002 for 
a discussion on the chemistry of the HCO molecule). 



Table 1. Parameters of gaussian fits of the SiO(2-l), H 13 CO + (1-0) 
and HCO(l-O) lines of Fig. 6. Errors (in brackets) are l<x values. 
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SiO(2-l) 


27 


(6) 


0.96 (0.15) 


18 


(5) 


34 (6) 




H 13 CO + (1-0) 


37 


(6) 


1.37 (0.16) 


17 


(5) 


35 (5) 
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SiO(2-l) 


22 


(5) 


0.54 (0.10) 


17 


(5) 


23 (5) 




H 13 CO + (1-0) 


43 


(5) 


1.15 (0.12) 


17 


(5) 


25 (5) 
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SiO(2-l) 


48 


(5) 


2.26 (0.21) 


43 


(5) 


45 (6) 




H 13 CO + (1-0) 


42 


(5) 


1.10 (0.14) 


42 


(5) 


25 (5) 




HCO(F=l-l) 


17 


(5) 


0.58 (0.16) 


50 
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32 (9) 




HCO(F=2-l) 


16 


(5) 


0.56 (0.16) 


49 


(6) 


32 (9) 
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SiO(2-l) 


31 


(5) 


1.61 (0.20) 


51 


(5) 


48 (8) 




H 13 CO + (1-0) 


20 


(5) 


0.52 (0.12) 


52 


(5) 


24 (5) 


E 


SiO(2-l) 


42 


(6) 


1.03 (0.11) 


13 


(5) 


23 (5) 




H 13 CO + (1-0) 


44 


(6) 


1.05 (0.12) 


15 


(5) 


22 (5) 



are, on average, a factor of 2 larger than the linewidths of 
H 13 CO + . The estimated velocity dispersion of the GMCs (<x v ) 
is ~ 20-25 km s _! , in the northern spiral arm, and ~ 10 km s , 
in the nuclear ring. This suggests that the degree of apparent 
turbulence measured on GMC-like scales in the SiO emitting 
gas is enhanced compared to that of the more quiescent dense 
gas component traced by H I3 CO + in this region. The latter im- 
plies that the SiO(2-l)-to-H I3 CO + (l-0) ratio measured at the 
wings of the SiO lines should be even larger than the velocity- 
integrated ratio Ri derived above for the northern spiral arm. 
This would only reinforce the case of an enhanced SiO chem- 
istry in gas phase in this region. 

Fig. 6 shows the SiO and H 13 CO + spectra observed towards 
the positions of GMCs A-to-E (the corresponding gaussian fits 
are listed in Tab. 1). As sketched in Fig. 6, we report on the 
tentative detection of two out of the four hyperfine line compo- 
nents of the Nk-k* - 1 -0 group of transitions of HCO towards 
the position of GMC C. These correspond to the rest frequen- 
cies 86.671 GHz (J = 3/2 - 1/2, F = 2 - 1) and 86.777 GHz 
(J = 1/2-1/2, F = 1-1). The HCO(F = 2-l)-to-H I3 CO + (l- 
0) intensity ratio inferred at GMC C is ~ 0.51 + 0.16; compa- 
rable ratios (~ 0.15 - 0.5) were derived by Garcia-Burillo et 
al. (2002) for the -650 pc nuclear disk of M 82, where the 
large HCO abundances LY(HCO)~ 4 x 10" 10 ) indicate that the 
whole inner disk can be viewed as a giant PDR. The tentative 
detection of HCO towards C would be well accounted if UV 
fields are partly driving the chemistry of the molecular clouds 
closest to the embedded star forming complex identified in the 



4. SiO fractional abundances 

We have estimated the column densities (N) of SiO and 
H 13 CO + from the PdBI map of IC 342 using a Large Velocity 
Gradient (LVG) code. Our aim is to estimate the abundance 
of SiO relative to H 2 (X(SiO)) inferred here from the SiO- 
to-H 13 CO + column density ratio (R N ). Values of TY(SiO) and 
7V(H 13 CO + ) are derived from measured integrated intensities, 
assuming a plausible range of physical conditions for the gas. 
Given that the emission of both species is optically thin, R?j is 
proportional to the SiO/H 13 CO + intensity ratio (#/). Since SiO 
and H 13 CO + have similar dipole moments and the observed 
transitions have comparable upper state energies, it is reason- 
able to assume the same physical conditions for both species. 
We have run five LVG models covering a range of gas densi- 
ties (n(H2)) from 10 4 cirr 3 to 10 6 cirr 3 . The explored interval 
encompasses the total range of molecular gas densities deter- 
mined in the GMCs of IC 342 from multitransition studies of 
CO and HCN (e.g., see Schulz et al. 2001). As gas kinetic tem- 
perature, we adopt a value — 50 K (Downes et al. 1992). 
The value of is not critical in the estimate of Rpj within the 
range of «(H2) explored in these calculations: the inferred col- 
umn density ratios are similar within the temperature interval 
20 K< T K <80 K and R N increases, at most, by 40% if T K is 
lowered to 10 K. 

We show in Fig. 7 the values estimated for R^ towards 
GMCs A-to-E for T K =50 K particularized for the different 
n(H2) values. On average, R^ is seen to increase by a factor 
of ~ 4 when «(H2) is lowered from 10 6 cirr 3 to 10 4 cm 4 . The 
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Fig. 4. Velocity channel maps of Si0(2-1) (grey scale) and H 13 CO + (1-0) (contour levels) in IC 342. Both scales are from 2.5<r to 10.5cr in 
steps of 1.5cr; \a = 1.1 mjy beanr' in a 10 km s~' channel. The channel LSR-velocities and the positions of GMCs A-to-E are indicated in 
each panel. 



fractional abundance of SiO is inferred from assuming a shows a fairly stable abundance in a large variety of phys- 

standard abundance for H 13 CO + . The assumption of a standard ical and chemical environments (see discussion in Garcfa- 

value for X(H 13 CO + ) is supported by observations of molec- Burillo et al. 2000). Here we adopt X(H 13 CO + )=2.5xl0" ll) ; 

ular clouds in our Galaxy. Contrary to SiO, for which mea- this corresponds to a typical abundance of the main iso- 

sured abundances are seen to differ by several orders of magni- tope X(H 12 CO + )=10~ 8 and to an isotopic ratio [ 12 C]/[ 13 C]^40 

tude between quiescent clouds and shocked regions, H 13 CO + (Henkel et al. 1998). 
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Fig. 5. SiO(2-l) and H 13 CO + (1-0) position-velocity (p-v) diagrams (b(top right) and c(bottom right)) taken along the ID-strips highlighted 
in a(left) (overlay of the SiO(2-l) and H 13 CO + (1-0) intensity maps as shown in Fig. 2b). The position angles of the p-v diagrams are chosen 
to maximize the contrast between the linewidths measured in the northern spiral arm and in the nuclear ring. In the two p-v plots, the SiO(2-l) 
brightness is represented in color scale (2.6 to 1 1 mJy beam -1 ) and white contours (from 3 to 1 1 in steps of 1.5 mJy beam -1 ), while H 13 CO + (l- 
0) brightness levels appear in black contours (same levels as above). The approximate location of GMCs B,C, D and E along the strips is 
indicated in b(top right) and c(bottom right). 
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Fig. 6. SiO(2-l) and H 13 CO + (1-0) spectra observed towards GMCs A-to-E. Spectra are smoothed to 10 km s -1 resolution. LSR-velocities are 
referred to the frequency of SiO(2-l). A horizontal arrow marks the expected velocity of H 13 CO + (1-0) at each panel according to the measured 
centroid for the SiO(2-l) line (i.e., redshifted — h-321 km s -1 from the SiO line). The expected velocities for the four hyperfine lines of HCO 
are marked by vertical arrows on the panel of GMC C. Vertical ±ltr errorbars indicate the expected peak temperatures of the HCO hyperfine 
lines, fixing the HCO(7 = 3/2 - 1/2; F = 2 — 1) temperature from observations and assuming optically thin emission for the four lines. 
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The estimated abundance of SiO is > 2 x 10" 10 for all 
GMCs, i.e., at least two orders of magnitude larger than the 
typical SiO abundances of Galactic quiescent clouds (Martin- 
Pintado, Bachiller & Fuente 1992). This lower limit on the 
beam-averaged value of X(SiO) measured here on GMC-like 
scales (~ 80 pc) indicates that shock chemistry is at work in 
the inner r ~ 320 pc disk of IC 342. Most remarkably, we 
see an overall N-S gradient in the value estimated for X(SiO). 
On average, the abundance of SiO is nearly one order of mag- 
nitude larger in the northern spiral arm (e.g., ~l-4xl0~ 9 at 
GMC D) than in the nuclear ring (e.g., ~2-7xl0~ 10 at GMC 
B) within the explored range of densities. Schulz et al. (2001) 
have estimated the average densities of GMCs A-to-E based 
on a multitransition study of HCN and their isotopes. The re- 
sults of this study indicate that average densities are a factor 2 
larger in the nuclear ring GMCs compared to that measured in 
the northern spiral arm. This would imply that the N-S gradi- 
ent on the value of X(SiO) may be even larger than estimated 
above. Moreover, the overall mass budget of IC 342 is seen 
to be heavily weighted by molecular gas with typical densities 
~ 10 3 - 10 4 cm" 3 (Downes et al. 1992; Israel & Baas 2003). 
Should this diffuse gas partly contribute to the emission of both 
species this would imply that the SiO abundances may have to 
be boosted for all the GMCs (see Fig. 7) compared to the values 
reported above, thus reinforcing the case for shock chemistry. 



n(H 2 )(cm 3 ) 




B A E C D 

Fig. 7. SiO-to-H 13 CO + column density ratios and SiO abundances in 
GMCs A-to-E derived from a LVG calculation, assuming r K =50 K 
and five different values of n(H2) from 10 4 to 10 s cirT 3 . GMCs are 
ordered along the X-axis by increasing X(SiO). 



The beam-averaged values of X(SiO) derived above are 
necessarily lower limits to the real abundance attained by SiO 
in the fraction of dense molecular gas which is being processed 
by shocks in IC 342 (hereafter X(SiO)U oc *). If we denote by 
f shock the a priori unknown fraction of shocked dense gas, 
X(SiO) can be formally factorized as: 

X(SiO) = f shock x X(SiO)| shock (1) 



The notable N-S gradient found in the derived value of 
X(SiO) between the nuclear ring and the northern spiral arm 
could be accounted by two extreme scenarios: first, an order of 
magnitude change in f s h oc k between the two regions, or alterna- 
tively, a similar change in X(SiO)| shock- While the first scenario 
would imply that the typical column densities of shocked gas 
in the northern spiral arm are ~ 6 - 7 times larger than in the 
nuclear ring, the second scenario would call for a notably dif- 
ferent shock velocity regime (v s hock) m the two regions. In this 
case v s hock should be significantly larger in the northern spiral 
arm. As is discussed in Sect. 7, finding the right scenario be- 
tween the two proposed above is key to shed light on the nature 
of the driving mechanism of shocks in IC 342. In Sect. 5 we 
compare the emission of SiO in IC 342 with that of CH3OH, 
another molecular shock tracer. In particular, we use the SiO- 
to-CH 3 OH ratio derived for IC 342 to explore the origin of 
shocks in this galaxy. 

5. Molecular shock chemistry in IC 342 

5.1 . Tracers of shocks in molecular gas 

The significant enhancement of SiO in gas phase is considered 
to be an indication that shock chemistry is at work in molecular 
gas (Martin-Pintado et al. 1992). The injection of Si-bearing 
material from dust grains into the gas phase, either through 
sputtering or grain-grain collisions, can explain the measured 
abundances of this molecule in shocked regions (Schilke et 
al. 1997; Caselli et al. 1997). Shocks are often invoked to ac- 
count for the large abundances of other molecular species mea- 
sured in bipolar outflows. This is the case of CH3OH (Bachiller 
et al. 1995). Although qualitatively similar, shocks character- 
ized by different velocity regimes are expected to process to 
a different extent dust grains in molecular gas. Fast shocks 
(vshock >15-20 km s~') can destroy the grain cores, liberat- 
ing refractory elements to the gas phase (Schilke et al. 1997; 
Caselli et al. 1997). In contrast, while slow shocks (v s hock <10- 
15 km s" 1 ) are not able to destroy the grain cores, they can 
heavily process the icy grain mantles. The assumed different 
location of Si-bearing material (cores) and solid-phase CH3OH 
(mantles) in dust grains makes of SiO and CH3OH good tracers 
of fast and slow shocks, respectively. Furthermore, for veloci- 
ties above ~ 10-15 km s _1 shocks could destroy the molecules 
in gas-phase of volatile species such as CH3OH (Garay et al. 
2000; J0rgensen et al. 2004). The dissociation of SiO by shocks 
would require velocities >50-60 km s _1 (i.e. J shocks), how- 
ever. 

If we consider both processes, i.e., the injection of grain 
material and the disruption of molecules in gas phase, we 
can conclude that an increase in the typical velocity regime 
of shocks (v shock) will certainly favour an enhancement of the 
abundance of SiO in the shocked gas (X(SiO)|.,/, OCJ t). This will 
be at the expense of increasing the SiO-to-CH30H abundance 
ratio in molecular gas. Therefore, a variation in the SiO-to- 
CH3OH intensity ratio can be taken as a evidence for a change 
in Vshock- A quantitative comparison of the SiO and CH3OH 
maps of IC 342 could thus help to discern if the typical shock 
velocity regime changes across the galaxy disk. 
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Fig. 8. We represent the SiO-to-CHsOH integrated intensity ratio as a 
function of the SiO-to-H 13 CO + ratio (R 7 ) in GMCs A-to-E. CH 3 OH 
data have been taken from Meier & Turner (2005). 



5.2. Tracers of shocks in IC 342: SiO and CH 3 OH 

Meier & Turner (2005) used the OVRO interferometer to map 
the emission of the 2k- It line of CH3OH in IC 342 with a 
resolution similar to that of the SiO map. This spatial resolu- 
tion is comparable to the typical GMC-scales (~ 80 pc). We 
have compared the emission of SiO(2-l) to that of CH30H(2k- 
lk) in order to derive the SiO-to-CHsOH intensity ratio in 
GMCs A-to-E. Fig. 8 represents these ratios as a function 
of Rj for the nuclear ring and northern spiral arm GMCs. As 
is shown in Fig. 8, the reported nearly order of magnitude 
change in Ri is not corresponded with a similar change in the 
I(SiO)/I(CH30H) ratio between the northern spiral arm and the 
nuclear ring. The I(SiO)/I(CH30H) ratio is fairly constant and 
close to ~0.14 1 for all GMCs, though we find tentative evi- 
dence for a larger value in GMC D (~0.24±0.07). Taken to- 
gether these results indicate that, at first order, the N-S gradient 
measured in /?/ between the northern spiral arm and the nuclear 
ring can be mostly attributed to a variation of f shock- 

The line ratios above are derived from velocity-integrated 
intensities and are beam-averaged on scales which are typi- 
cal of GMC-like units at the distance of IC 342. However, the 
analysis of the line profiles of SiO, CH3OH and H 13 CO + can 
provide information on the shock velocity regime on scales 
smaller than the beam. Fig. 9 displays the SiO-to-H 13 CO + 
and SiO-to-CH 3 OH velocity-width ratios derived for GMCs 
A-to-E. These ratios show a different behaviour in the nu- 
clear ring and in the northern spiral arm. In the ring, linewidths 
and velocity centroids for SiO, H 13 CO + and CH3OH are vir- 
tually identical. In contrast, SiO lines are a factor of ~2 wider 
than that of H 13 CO + in the northern spiral arm, as reported in 
Sect. 3.2.2. CH3OH lines represent a case intermediate between 
these two extremes: SiO lines are a factor of ~1.4 wider than 



1 Assuming the CH 3 OH column densities given by Meier et 
al. (2005) these intensity ratios translate into SiO-CHsOH column 



density ratios of ~ 0.01 
larger if n ~ 10 4 cm -3 ). 
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Fig. 9. SiO-to-H 13 CO + and SiO-to-CH 3 OH velocity width ratios 
measured for GMCs A-to-E. These ratios are represented as a func- 
tion of the SiO-to-H 13 CO + integrated intensity ratio, Rj. 



that of CH3OH in the spiral arm. The differences between SiO, 
CH3OH and H 13 CO + are evident in the linewidths, but veloc- 
ity centroids are the same within the errors. As discussed in 
Sect. 7, this suggests that the apparent turbulence of shocked 
molecular gas is enhanced compared to the more quiescent gas. 

While it seems that the enhancement of f shock in the north- 
ern spiral arm explains the bulk of the reported increase of 
X(SiO) in this region, the observed differences in the line pro- 
files of SiO and CH3OH suggest that a fraction of the shocked 
gas in the arm presents a higher X(SiO)\ s hock- This implies that 
v shock would be a factor of ~2 larger in the spiral arm region. 

6. SiO emission in Galaxies 

There is ample observational evidence that SiO thermal emis- 
sion can be locally enhanced in star forming molecular clouds 
of our Galaxy (Martin-Pintado et al. 1992; Bachiller et al. 
2001). More recently, SiO has been revealed to be a tracer 
of shock chemistry also in galaxy nuclei, including our own 
Galaxy (Sage & Ziurys 1995; Martin-Pintado et al. 1997; 
Garcfa-Burillo et al. 2000, 2001b; Usero et al. 2004). We re- 
vise below the various driving mechanisms which have been 
proposed to explain the onset of shock chemistry in the Galaxy 
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and in galaxies in general, and study their applicability to the 
nucleus ofIC 342. 

6.1. SiO emission in the disk of our Galaxy 

SiO thermal emission is observed towards star-forming clouds 
in the disk of our Galaxy. In particular, the strongest SiO emit- 
ters are the bipolar outflows located around protostellar objects 
(Martin-Pintado et al. 1992; Bachiller et al. 2001). In the first 
stages of the star formation process, bipolar outflows interact 
with the ambient molecular gas, inducing molecular shocks 
which are able to increase the abundances of some molecular 
species (like SiO and CH3OH) by several orders of magnitude 
with respect to quiescent gas (Bachiller et al. 1995, 1997). 

The abundance of SiO in IC 342 reaches the largest value in 
the northern spiral arm (Sect. 4), where the evidence of ongoing 
star formation is scarce (Sect. 3.1). However, it could be argued 
that SiO emission in IC 342 is probing the deeply embedded 
phase of a young star formation episode. The episode would 
be spread on scales of ~a few 100 pc and would not visible in 
Ha or thermal radio-continuum. Although with these tight con- 
straints this explanation is unlikely, we can discard it on more 
quantitative grounds comparing the SiO emission in IC 342 
with a subset of bipolar outflows of the Galaxy (Tab. 2). For the 
reasons explained below, we have purposely selected outflows 
for which there are maps of both SiO(2-l) and CH 3 OH(2 k -l k ) 
available in the literature. Furthermore, since SiO luminosities 
are larger for more massive and younger objects, the outflows 
have been chosen to cover a wide range in mass and age of the 
protostars. It is assumed that the largest SiO-to-CHjOH lumi- 
nosity ratios correspond to the less evolved objects (Bergin et 
al.1998). 

The mean SiO intensity of the northern spiral arm and the 
nuclear ring of IC 342, derived inside the region defined by 
the lowest contour of Fig. 2, is ~0.93 K km s _1 . We have then 
inferred the surface density of outflows which would be re- 
quired to reproduce the observed SiO intensity in IC 342. As 
is shown in Tab. 2, this number density of outflows range from 
~0.4 to ~12 outflows pc~ 2 . Though in star-forming regions of 
our Galaxy like OMC2/3 or NGC 2068, surface densities of 
a few outflows pc~ 2 have been reported on scales of 1-2 pc 2 
(Reipurth et al. 1999; Mitchell et al. 2001), it is very unlikely 
that these average surface densities can be attained in the SiO 
disk of IC 342, which is several hundred pc 2 in extent. 

Furthermore the star forming rate (SFR) of IC 342 poses 
tight constraints on the upper limit to the expected density of 
outflows. The SFR in a galaxy can be easily estimated from 
the FIR luminosity (Kennicutt 1998). In the case of IC 342, 
L FIR = 1.5 X 10 9 L.0 in the central r ~ 15" (Becklin et al. 
1980; scaled to D=3.3 Mpc). Assuming that the northern spi- 
ral arm and nuclear ring are the main contributors to the 
SFR of IC 342, we estimate for this region a SFR density of 
3 x 10~ 6 M0 yr~' pc~ 2 . From the SFR density above, we can 
then derive an upper limit to the surface density of outflows, as- 
suming a Salpeter law for the IMF (d log(A0/d log(m) = -2.35 
over m = 0.1-100 Mq) and a timescale for the pre-stellar 
phase of ~ 10 4 yr. The derived upper limit to the total density 



of outflows provided by the SFR in IC 342 is ~ 9 x 10~ 2 out- 
flows pc 2 . Even for outflows like CB3 (or any other SiO lumi- 
nous YSO), the required density would be ~5 times larger than 
that provided by the SFR of IC 342. Moreover, massive young 
bipolar outflows like CB3 should be the minority among YSOs: 
for a standard Salpeter IMF, only 4% of the objects would have 
masses above 1 Mq. In addition, the observed SiO-to-CHsOH 
average ratio of IC 342 (~ 0.14; see Fig. 8) is similar to YSOs 
which are much more evolved than CB3. 

In summary, we can discard the interpretation of the large- 
scale shocks in IC 342 in terms of a collection of YSOs as- 
sociated with an embedded star formation episode in the inner 
r ~ 320 pc disk of the galaxy. 

6.2. SiO emission in galactic nuclei 
6.2.1 . The center of our Galaxy 

The nucleus of our Galaxy shows widespread SiO emission 
that, in contrast to that observed in the Galactic Disk, is 
not related to recent star formation. The first large-scale SiO 
observations of the Galactic Center (GC) of Martin-Pintado 
et al. (1997) detected the emission of the SiO(l-O) line in 
a ~ 150 pc-diameter circumnuclear disk. In a later paper, 
Hiittemeister et al. (1998) detected the SiO(2-l) emission in 32 
GC molecular clouds located inside a ~300 pc-diameter disk 
which extends from Sgr C to a position at slightly higher posi- 
tive longitudes than Sgr B2. High SiO abundances are derived 
for these clouds (~ several 10~ 10 - several 10~ 9 ), indicative 
of shock chemistry. This scenario has received further support 
from recent observations revealing the complex alcohol chem- 
istry of GC clouds (Martin-Pintado et al. 2001; Requena-Torres 
et al. 2005, in prep.). In these clouds the large abundances of 
ethanol (C2H5OH) and CH3OH evidence the erosion of dust 
grain mantles on large scales. 

Besides the long reported similarities between our Galaxy 
and IC 342 (e.g., Downes et al. 1992), the inner few hun- 
dred pc in the two galaxies seem to be the scenario of 
large-scale molecular shocks leading to dust grain processing. 
Furthermore, the efficiency of shocks appear to be comparable 
in both galaxies. Despite the different spatial resolution of SiO 
observations in IC 342 (~80 pc) and in our Galaxy (~2 pc), 
the derived SiO abundances are similar in the two objects. 
Moreover, the SiO-to-CHsOH abundance ratios are also similar 
in the GC clouds (~ 0.01 - 0.03) 2 and in IC 342 (~ 0.01 - 0.04; 
see Sect. 5.2). 

These similarities taken together, we can hypothesize that 
the mechanism explaining the onset of large-scale molecular 
shocks in the GC and in IC 342 is likely to be the same. Shocks 
identified in GC clouds, also unrelated to ongoing star forma- 
tion, have been attributed different causes, however: the inter- 
action with GC non-thermal filaments, with supernovae rem- 
nants or with the expansion bubbles created by Wolf-Rayet 
stars have been discussed by Martin-Pintado et al. (1997). 

2 Calculated from SiO and C 2 H 5 OH column densities measured in 
GC clouds (Martin-Pintado et al. 1997, 2001), assuming a constant 
C 2 H 5 OH-to-CH 3 OH abundance ratio of ~0.05 (Requena-Torres et 
al. 2005, in prep.). 
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Table 2. SiO(2-l) and CH 3 OH(2 k -l k ) luminosities of a sample of bipolar outflows. Col. 2: mass of the class object. Col. 3: luminosity of 
the class object. Col. 4: SiO luminosity integrated within the SiO emitting region of the outflow. Col. 5: SiO-to-CH 3 OH luminosity ratio. 
CoJ. 6: surface density of outflows required to obtain the mean /(SiO) in the northern spiral arm+nuclear ring of IC 342. Col. 7: References: a: 
Bachiller et al. (2001); b: Gueth et al. (1997); c: Fuente et al. (2005b); d: Codella & Bachiller et al. (1999). 
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Hiittemeister et al. (1998), using observations that extended 
over a wider region in the GC, found the largest SiO abun- 
dances where the likelihood of cloud-cloud collisions, induced 
by the Galactic bar potential, is the highest. As discussed in 
Sects. 6.2.2 and 7, the latter scenario is the preferred one in 
IC 342. 

6.2.2. The nuclei of external galaxies 

Mauersberger & Henkel (1991) detected the emission of 
SiO(2-l) in the starburst galaxy NGC 253, using the 30 m 
IRAM telescope with a resolution of ~28". This detection, the 
first reported for SiO outside the Milky Way, was followed by 
a 9-galaxy survey made with the NRAO 12 m antenna by Sage 
& Ziurys (1995) with ~67" spatial resolution (equivalent to 
~0.5-2 kpc). In this survey, where 5 galaxies were detected 
in SiO, Sage & Ziurys (1995) found no correlation between the 
abundance of SiO, characterized in their work by the SiO/N2H + 
ratio, and the efficiency of star formation (given by the SFR 
per unit dynamical mass, i.e., SFR/M</ V „). This surprising result 
was at odds with the classical framework where shock chem- 
istry is driven by YSOs in star forming regions of the Galactic 
disk. The global SiO abundance measured on 0.5-1 kpc-scales 
is seen to vary significantly among starburst galaxies: it can 
reach ~10~ 9 in IC 342, whereas it is 1/20 of this value in M 82 
(see Table 3 and Martin-Pintado et al. 2005, in prep.). 

The occurrence of large-scale molecular shocks may arise 
at different stages during the typical lifetime of a star- 
burst episode (Rieke et al. 1988; Garcia-Burillo & Martin- 
Pintado 2001a). In the pre-starburst phase (I), density wave in- 
stabilities induce gravitational torques and drive the infall of 
large amounts of gas towards the nucleus. Large-scale shocks 
may be at work related to an enhanced compression of gas 
and an increased rate of cloud-cloud collisions in the potential 
wells of spiral arms and/or bars. Once the first massive stars are 
formed in a second phase (II), bipolar outflows can produce lo- 
cally molecular shocks in YSOs. In a later stage, corresponding 
to an evolved starburst, the elevated rate of SN explosions may 
lead to the disruption of the disk during the expansion of the 
so-called hot bubble. Episodes of mass injection from the disk 
into the halo could be accompanied by molecular shocks. 



The advent of high-resolution SiO images has been key to 
help discern the different sources of shock chemistry in ex- 
ternal galaxies and thus identify an evolutionary path along 
the starburst sequence depicted above (see Table 3 and refer- 
ences therein). In particular, the puzzling result issued from 
first single-dish SiO surveys starts to be understood when ob- 
servations allow us to zoom in on molecular galaxy disks on 
scales <100 pc: 

- Large-scale shocks near the outer Inner Lindbland 
Resonance (oILR) of the NGC 253 stellar bar can account 
for the outer SiO disk of this galaxy, which extends well 
beyond the nuclear starburst. The outer disk SiO emission 
gives away the pre-starburst phase (I) in NGC 253 (Garcfa- 
Burillo et al. 2000). More recently, the detection of S0 2 , NS 
and NO emission and the analysis of the sulfur chemistry 
in NGC 253 have confirmed that shocks are at work in the 
nucleus of this galaxy (Martin et al. 2003, 2005). 

- Outflows driven by YSOs during phase II can partly ex- 
plain the measured average abundances of SiO in the inner 
disk of NGC 253 (Garcia-Burillo et al. 2000). However, the 
bulk of the SiO emission likely stems from molecular cloud 
shocks in the inner Inner Lindbland Resonance (ilLR) of 
the NGC 253 stellar bar. 

- The emission of SiO extends noticeably out of the galaxy 
plane in M 82, tracing the disk-halo interface where 
episodes of mass injection from the disk are building up 
the gaseous halo (phase III) (Garcia-Burillo et al. 2001b). 
The PdBI maps of M 82 made in SiO and HCO illustrate 
how two different gas chemistry scenarios can be simulta- 
neously at play in the same galaxy though at different loca- 
tions: shocks in the disk-halo interface and PDR chemistry 
in the galaxy disk which hosts an evolved starburst (Garcia- 
Burillo et al. 2002; Fuente et al. 2005a). 

The case of the Seyfert 2 galaxy NGC 1068 has been stud- 
ied by Usero et al. (2004) who estimate a large abundance of 
SiO in the r ~ 200 pc CND of the galaxy (X(SiO)~a few 10~ 9 ). 
The enhancement of SiO is attributed by Usero et al. (2004) to 
the evaporation of very small (10 A) silicate grains (VSG) by 
X-rays. Alternatively, silicon chemistry could be also driven by 
pre-starburst shocks related with the density wave-resonances 
of the CND (Garcia-Burillo et al. 2005, in prep.) 
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Table 3. SiO emission in external galaxies. Col. 2: SiO abundance averaged within 28" (from observations with the IRAM 30 m telescope); 
the averaging scale is ~ 0.5 kpc, for NGC 253, IC 342 and M 82, and ~ 2 kpc for NGC 1068. Col. 3: components resolved with the PdBI 
contributing to the global SiO emission; CoJ. 4: SiO abundance in the resolved components. Col. 5: mechanism proposed to be enhancing the 
SiO abundance; underlying physical process, in brackets (DW stands for 'density waves'). Col. 6: References: a: Garcfa-Burillo et al. (2000); 
b: this work; c: Garcfa-Burillo et al. (2001b); d: Usero et al. (2004); e: Martin-Pintado et al. (2005, in prep.). 
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As extensively argued in Sect. 6.2.1, the SiO emission in 
IC 342 cannot be explained by ongoing star forming activity 
(phase II). The M 82 scenario (phase III) can be also ruled out 
in IC 342. SiO emission in the disk of IC 342 extends well be- 
yond the distribution of supernovae remnants (SNRs) (Condon 
et al. 1982; Bregman et al. 1993). Furthermore, the low super- 
novae rate of IC 342 ( < 0.04 yr _ 1 for D=3.3 Mpc; Condon et 
al. 1982) yields a ~ 70 times smaller energy deposition by the 
SN of IC 342 compared to M 82. Finally, X-rays are not ex- 
pected to be a dominating agent in the chemistry of molecular 
gas in IC 342 either. Nearly ~35% of the observed hard X-ray 
emission in IC 342 (2 keV< E <10 keV) comes from a circum- 
nuclear disk of r ~ 8" (Bauer, Brandt & Lehmer 2003). The 
hard X-ray luminosity of this source is three orders of magni- 
tude lower than that of NGC 1068 (Ogle et al. 2003), however. 
Furthermore, the 6.4 keV Fe Ka line, which probes the pro- 
cessing of neutral gas by X-rays, is absent from the spectrum 
ofIC342. 

The exclusion of all the alternative explanations leads us 
to conclude that the large-scale shocks identified in the in- 
ner r ~ 320 pc of IC 342 arise in the pre-starburst phase (I) 
(see Sect. 8). The following Section discusses the efficiency of 
density-waves at producing shocks in molecular gas in IC 342. 

7. Density waves and shocks in IC 342 

The IC 342 bar shapes the distribution and kinematics of 
molecular gas in the central r ~ 320 pc of the galaxy (Turner 
& Hurt 1992; Schinnerer et al. 2003; Meier & Turner 2005; 
this work). The spiral-like morphology of the molecular disk 
and the detection of non-circular motions (> 50 - 60 km s" 1 in 
the northern spiral arm, deprojected onto the galaxy plane) are 
reminiscent of the typical bar-driven dynamics. The SiO abun- 
dances measured in the inner molecular disk of IC 342 proves 



unambiguously that the bar is producing large-scale molecular 
shocks. However, the detection of SiO emission constrains the 
velocity regime of shocks to lie between > 15-20 km s -1 (for 
grain cores to be significantly disrupted) and < 50 - 60 km s -1 
(to prevent dissociation of SiO molecules), i.e., a velocity range 
which lies significantly below the lower limit set to the non- 
circular motions measured across the IC 342 bar. Furthermore, 
the detection of CH3OH emission across the bar suggests that 
the shocked molecular gas emitting in CH3OH cannot be char- 
acterized by v shock > 50 - 60 km s -1 . These observational con- 
straints imply that the input kinetic energy provided by stream- 
ing motions must be first dissipated making room for a lower 
velocity regime that corresponds to the emission of the molec- 
ular shock tracers observed across the IC 342 bar. 

Large-scale shocks driven by density waves have been 
long predicted by numerical simulations of spiral/barred galax- 
ies, made following either hydrodynamical schemes (Roberts 
1969; Athanassoula 1992) or ballistic ones (Casoli & Combes 
1982; Combes & Gerin 1985). In ballistic models, which likely 
provide a more realistic representation of the clumpy dense 
ISM in galaxies, molecular shocks should arise subsequently 
after cloud-cloud collisions. The number of collision events is 
enhanced by an increase of orbit crowding along the potential 
well of the bar. As noted in Sects. 3.2.2 and 5.2, SiO, H 13 CO + 
and CH3OH lines have all similar velocity centroids, but dif- 
ferent linewidths over the spiral arm region. The largest widths 
correspond to SiO (whose lines are a factor of 2 larger than that 
of H 13 CO + ), with CH3OH representing an intermediate case. 
This suggests that the apparent turbulence of shocked molec- 
ular gas, traced by SiO and CH 3 OH, is enhanced compared to 
the state of the more quiescent dense gas medium (traced by 
H 13 CO + ). This effect, especially relevant in the northern spi- 
ral arm region (see Fig. 5), suggests that the molecular shocks 
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arise at a stage of turbulent dissipation and not during the early 
phase of the encounter, when cloud-cloud relative velocities 
may lead to the dissociation of SiO and CH 3 OH. 

Assuming that the internal structure of the colliding molec- 
ular clouds is highly clumpy (e.g., Falgarone & Puget 1985), 
it is plausible to assume that a fraction of the kinetic energy 
dissipated during a cloud-cloud collision can cascade down to 
smaller scales. This would increase the turbulent motions of 
clumps composing the end product of any cloud-cloud colli- 
sion. On the simulations front, Kimura & Tosa (1996) have in- 
deed found indications that the internal turbulence of clumpy 
colliding molecular clouds can increase after an encounter (see 
also Bonnell et al. 2005). The input kinetic energy typically 
involved in a cloud-cloud collision should be much larger in 
the northern spiral arm than in the nuclear ring. Once dissi- 
pated, this different input energies would end up producing a 
higher turbulence in the shocked molecular gas of the spiral 
arm (cr v ~ 20- 25 km s _1 ; see Sect. 3.2.2) compared to the 
nuclear ring (cr v ~ 10 km s _1 ; see Sect. 3.2.2). As observed, 
in this scenario the largest difference between these two re- 
gions would be the amount of shocked molecular gas mass 

(fZck ~ 5 - 1X fTock> See SeCt 5 - 2 > and n0t the Sh0ck Vel0C " 

ity regime which is here equal to cr v (cr^™ ~ 2 x cr™ g ; see 
Sect. 5.2). 

The fraction of molecular gas actually involved in the pre- 
starburst shocks produced by the IC 342 bar defines the rele- 
vance of this process. We can estimate this fraction in IC 342 
making use of Eq. 1 . Given the similar properties of shocks in 
the GC and IC 342 ( Sect. 6.2.1), we can reasonably adopt a 
value of X(SiO)\ s hock in IC 342 similar to that estimated in the 
GC clouds on ~l-2 pc-scales by Hiittemeister et al. (1998): 
X(SiO)L oc * > 10~ 8 . From the values derived for X(SiO) in 
IC 342 on ~ 80 pc scales, we conclude that f s hock ^ 0.1 
in the northern spiral arm and < 0.02 in the nuclear ring. 
In Appendix A, we estimate the rate of energy dissipated by 
shocks in the gas over the spiral arm region of this galaxy. 
While this estimate is rather approximate, the outcoming pic- 
ture underlines the potential role of large-scale shocks at drain- 
ing energy from the gas inflowing towards galactic nuclei. 

8. Conclusions and perspectives 

The high-resolution images showing the emission of SiO in the 
inner r ~ 320 pc disk of IC 342 reveal the onset of large-scale 
molecular shocks driven by the bar potential of this galaxy. The 
variation of the estimated SiO abundance inside the mapped 
region (from ~10~ 10 to ~10~ 9 ) and the comparison with other 
molecular tracers (CH 3 OH and H 13 CO + ) indicate that shocks 
process with uneven efficiency (see Sects. 4 and 5) the molec- 
ular gas reservoir of IC 342. The shocks seem to arise dur- 
ing cloud-cloud collisions at the stage when kinetic energy has 
partly dissipated in turbulent motions. The mass of molecular 
gas locally involved in shocks over the spiral arm region of 
IC 342 could amount to 10% of the dense gas reservoir. Taken 
together, these results underline the relevant role that large- 
scale molecular shocks can play at shaping the evolution of gas 
disks. 



These observations illustrate the occurrence of molecular 
shocks in galaxies which are not related with ongoing star for- 
mation. SiO emission in IC 342 traces a pre-starburst phase in 
molecular gas. Pre-starburst shocks have also been identified 
to be responsible of the intense ro-vibrational H2 lines detected 
at 2 ^m in 2 prototypical mergers: Arp 220 and NGC 6240 
(Rieke et al. 1985). More recently, Haas et al. (2005) have re- 
ported the detection of widespread emission of the v=0-0 S(3) 
line of H2 at 9.66 //m in the overlap region of the Antennae 
galaxy pair. Haas et al. (2005) interpret this emission as a tracer 
of shocks that will give rise to the first generation of stars in 
the region. In the case of the Antennae, Haas et al. (2005) hy- 
pothesize that shocks may not be the result of direct collisions 
of molecular clouds but arise instead from an overpressured 
medium that remains of the collisions between Hi clouds (Jog 
& Solomon 1992). 

In contrast to the Antennae, the bulk of the shocked gas in 
IC 342 is not expected to produce on-site star formation due to 
the inhibiting action of strong shear over the spiral arms. SiO 
emission in IC 342 traces the sites where molecular clouds dis- 
sipate a fraction of their energy through collisions; the energy 
lost helps the gas to fall to the nuclear ring where it will feed 
the starburst. The SiO map of IC 342 provides a snapshot view 
of the pre-starburst phase during the fueling process driven by 
density waves. Higher-resolution observations are required to 
provide new constraints on the details of how density waves 
operate to produce molecular shocks. In particular, we expect 
that the efficiency of shocks changes transversally to the spi- 
ral arms, as cloud orbits are re-oriented by the bar potential. 
Because of its closeness, favourable orientation and well de- 
fined spiral pattern, IC 342 is a good target for follow-up stud- 
ies. 

High-resolution SiO imaging is key to discern the different 
sources of shock chemistry which are activated at different lo- 
cations and at different moments in galaxy disks during a star- 
burst event. Being more than a mere tracer of exotic chemistry, 
SiO allows to probe unambiguously the regions where dust 
grains are being destroyed in galaxies due to the action of den- 
sity waves, star formation, galactic outflows or X-rays (Garcia- 
Burillo et al. 2000, 2001b; Usero et al. 2004). The study of the 
feedback influence of these phenomena in nearby galactic disks 
is paramount to constrain models of evolution and formation of 
galaxies at higher redshifts. 
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Appendix A: Energy dissipation by large-scale 
molecular shocks in IC 342 

Stellar bars can remove energy and angular momentum from 
the gas through gravitational torques and shocks. Gravity 
torques from the stellar potential on the gas are expected to 
be the most relevant drivers of gas inflow in galaxies. Torques 
created by large-scale stellar bars alone or helped in due time 
by other secondary mechanisms, such as nested stellar bars, dy- 
namical friction or viscosity, can drive gas inflow and feed the 
central activity of the inner 1 kpc playground in galactic nu- 
clei (e.g., Garcfa-Burillo et al. 2005). Large-scale shocks may 
also contribute significantly to the loss of energy. In this ap- 
pendix, we make a rough estimate of the rate of energy losses 
produced locally by the onset of large-scale shocks in the gas 
over the spiral arm region of IC 342. 

A shock-front propagating in a gas medium dissipates 
ordered kinetic energy into heat. For large Mach numbers 
(v shock » sound speed), the kinetic energy lost per unit 
shocked mass is ~ v^ hock /2, both in radiative and non-radiative 
shocks, (e.g. Draine & McKee 1993). The rate of energy dissi- 
pated per unit total mass due to shocks would be: 



/ 



shock 



TSiO 



(A.l) 
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/ is the total fraction of shocked gas (including the dense 
molecular component traced by SiO and the more diffuse 
gas), while T Si0 is the depletion timescale of SiO onto grains 
(~ 10 4 yr; Bergin et al. 1998). Since gas-phase SiO depletes 
quickly onto grains, the shocked gas traced by SiO must have 
been shocked within a time rsio- The expected values of / 
range between f shock (if dense and diffuse gas are assumed to 
be shocked to the same extent) and ~ 0.1 x f shock (if we as- 
sume that only the dense gas is preferentially shocked; in this 
case we take a fraction of dense gas of ~ 10% in IC 342 from 
Schulz et al. 2001). Assuming a typical v shock of 25 km s _1 , 
e s = -(0.3 - 3) x 10 3 km 2 s~ 2 Myr -1 for a corresponding range 
in/ = 0.01 -0.1. 

Assuming quasi-circular motions, we can estimate a dissi- 
pation timescale, t s , for the shocks to drain the specific energy 
of the gas, e (e v 2 /2 + (/>, where is the azimuthal velocity 
and (f>{r) is the mean gravitational potential): 

t s = -- x ^- (A.2) 

&s "arm 

The factor (2n)/(9 arm ) (where 9 arm is the angular width of 
the spiral arms at a given radius) accounts for the fact that, 
along an orbit, shocks dissipate energy only during the time 
spent by the gas in the spiral arm. Using 13 CO data, Turner & 
Hurt (1992) fitted a rotation curve in the inner 320 pc of IC 342 
of the form v (c (r) = a( Vl + br - 1) (with a, b constants); from 
the equation of motion, v 2 /r = d<p(r)/dr, we can estimate (f>(r) 
and e. The value of e results to be ^ (2 - 3) x v 2 /2. At the radius 
of GMC D (r ~ 220 pc) V(f ~ 60 km s" 1 and 6 arm l{2n) ~ 0.1, 
so we obtain e ^ 5 x 10 3 km 2 s~ 2 . The corresponding t s is 
~ 16-160 Myr, i.e. -0.7-7 rotations at the location of GMC 
D. 

Our estimates suggest that the energy could be drained effi- 
ciently from the gas by large-scale shocks along the spiral arms. 
While it is true that the inflow of gas is mostly constrained by 
the angular momentum transfer, rather than by the energy dis- 
sipation rate (draining angular momentum is more difficult), 
large-scale shocks could have indeed a non-negligible influ- 
ence in the dissipation of energy of the gas on its way to the 
nucleus. 



